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Abstract 

In a previous paper I showed how the ideal SLAG derivative and second-derivative 
■ operators for an infinite lattice can be obtained in simple closed form in position space, 

CN . and implemented very efficiently in a stochastic fashion for practical calculations on 

finite lattices. In this second paper I show how the small (order 1/N) errors introduced 
^ ' by truncating the operators to a finite lattice may be removed by a small adjustment 

•/^ ■ of coefficients, without incurring any additional computational cost. The derivation of 

i these results is again presented in a simple, pedagogical fashion. 

o 

O I 1. Truncating the SLAC derivative operators to fit on a finite lattice 

^ ' In a previous paper [1] we looked at the ideal "SLAC" specification of Drell, Weinstein 



and Yankielowicz [2] for the spatial derivative operator on an infinite one- dimensional lattice, 



QJ . which demands that —i times the Fourier transform of the derivative operator take on its 



ideal functional form, 

Pideal = P, 



X 

H 

\ within the first Brillouin zone. Such an operator avoids the pathologies of the fermion 
doubling problem or the violation of chiral invariance that have plagued other sub-optimal 
definitions of the derivative operator, but has a representation in position space that is very 
"nonlocal", 

1 r 1 . . 1 „ . 1 



^ideai/(x) = -{••■- -:fix + 4a) + -fix + 3a) - -fix + 2a) 
a L 4 6 2 



+ fix + a) - fix -a) (1) 
+ ^fix - 2a) - ^fix - 3a) + ^/(x - 4a) - ... }, 

and indeed links the lattice site at which we wish to take the derivative to every other site 
in the (one-dimensional) lattice. This has, historically, presented a barrier to the practical 
implementation of the SLAC operator in any large-scale lattice calculations. 

In [1] I proposed a solution to this computational barrier, namely, that the SLAC deriva- 
tive operator can be implemented in a stochastic fashion, so that for a one-dimensional finite 
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lattice of N sites the average number of computations required to implement each derivative 
operation is of order log N rather than of order N. 

However, even though the use of the ideal SLAC specification ensures that the repre- 
sentation (1) of the derivative operator is the best that can possibly be done for a given 
finite lattice spacing distance a, the "truncation" of the operation to "fit" on a lattice with 
a finite number N of sites in [1] introduces small errors (of order that detract from 

the "optimality" of the operator. It is desirable to remove these errors, to ensure that the 
stochastic SLAC derivative operator for a finite lattice is "optimal" — namely, the best that 
can possible be done. 

It is first worth seeing explicitly how the "truncation" operation detracts from the fidelity 
of the ideal SLAC derivative operator. Let us take a truly extreme example — a lattice of 
only one site! — and let us assume that the underlying continuum formalism specifies that 
we take the second-derivative of some function f{x). A simple derivation in [1] showed 
how one can obtain from first principles the expression first found by Drell, Weinstcin and 
Yankielowicz [2] for the ideal SLAC second-derivative operation in position space, for the 
case of an infinite lattice: 

^!dli/(^) = S{ ■ ■ ■ - + 4«) + Ifi^ + 3«) - \f{^ + 2a) 
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+ f{x + a)--f{x) + f{x-a) (2) 

- \f{x - 2a) + \f{x - 3a) - ^f{x - 4a) + ... |. 

What happens when we "truncate" this operation to apply to the extreme case of a lattice 
with just one site? Clearly, we will be left with only the middle term, namely. 



^ truncated /(-^y 



= -^/(O). (3) 
x=o 3a^ 



Assuming /(O) ^ 0, we have somehow estimated the second-derivative of a function evaluated 
at just one point in space to be some nonzero value! Of course, 1 said above that the errors 
in truncating the operator are of order 1/N; and here = 1, so the errors are of order unity; 
in other words, the value itself has absolutely no rehabihty at all. But is this the best we 
can do? 

I believe that it is not. When we restrict any physical formalism to a finite "box", 
we generally assume that the physical functions of relevance shall have periodic boundary 
conditions applied to them. (Indeed, it would be difficult to formulate a tractable analysis of 
most systems in momentum space if this stipulation were not to be made.) This is equivalent 
to thinking of a system of infinite spatial extent, made up of "boxes" stacked side-by-side 
along the x-axis, with the extra proviso that the contents of each box are to be considered to 
be identically equivalent to that of the "central box". If this is the case, then in the extreme 
case of a lattice with only one site, we are effectively stipulating that the value of any function 
at every site of the imagined infinite lattice (one such site in each of these identical "boxes" ) 
must be identical to the value of the function in the "central box" (centred on x = 0) . In 
other words, the function is effectively a constant on all of the (infinite number of) lattice 
sites: 

f{xn = no) — /(O) = constant. 
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In this context, our obtaining a nonzero estimate of the second-derivative in (3) looks very 
poor: surely, the second-derivative of a constant function is just zero! 



2. Applying periodic boundary conditions to the SLAC derivative operators 

The reason that the result (3) obtained above is so poor is simply because we have (by 
assumption) applied periodic boundary conditions to the function defined on the single lattice 
site, but we have not applied periodic boundary conditions to the SLAC derivative operator; 
rather, we simply "truncated" it when we had gone once around the lattice (in this case, after 
evaluating it at just the single lattice point). Now, the SLAC second-derivative operator (2) 
links the site in question to every other site on an infinite lattice. Applying periodic boundary 
conditions to it implies that we must go around and around the lattice an infinite number 
of times, calculating all of these infinite number of terms that contribute to the sum. 

Performing an infinite number of calculations does not sound very palatable, from a com- 
putational point of view. However, by the very assumption of periodic boundary conditions, 
the function f{x) itself is identically the same every time we "go past" it, and so it effectively 
factorises out of the sum, namely. 



= co/(0), 

x=0 



where 



2 f 1 1 1 7r2 ^ 1 1 1 1 

^°"^V-- 16 + 9-4 + ^- 6 +^-4 + 9-T6 + -t 

But this is simply an infinite sum of numerical factors, which needs to be only done once. 
Specifically, (4) is equivalent to 

Fortunately, the mathematicians have worked out this infinite sum for us, namely, 

^2 - "ly^) - 12 



E^-.(2)-;,. (5) 



n=l 

Thus, we find 

Co = 0, 

and hence the ideal SLAC second-derivative of any function defined on a lattice with just 
one site, with periodic boundary conditions applied to both the function and the operator, 
is identically zero, as argued above. 



3. A more realistic example 

Let us now consider a more realistic example. For simplicity, let's consider a one- 
dimensional lattice with, say, N — 100 sites. Let us also go back to considering the first- 
derivative of some function f{x) defined on the lattice. Now, since we are assuming periodic 
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boundary conditions, it doesn't really matter how we label the 100 lattice sites relative to 
the site that we wish to compute the derivative at; but, for definiteness, let us label them 
from n = —49 through n — +50 inclusive, with n — being the site that we wish to compute 
the first-derivative at. Now, the SLAC derivative operator (1) tells us to compute differences 
for increasing distances from n = 0, weighting the difference at distance n by {—l)"''^^/na. 
In [1], we simply assumed that we should stop when we get to the difference at n = ±49, 
since there is no n = —50. Let us consider, as an arbitrary example, the lattice position at 
n — +43. The contribution to the sum (1) from this position, according to the truncated 
prescription of [1], is simply 



43a 



We now wish to apply periodic boundary conditions to this derivative operation. Let us 
continue to concentrate on the position at n — +43. When the derivative operation "runs" 
through the lattice sites a second time, we will return to this same position, n = +43, after 
going a further distance oi N = 100 lattice sites. This implies a contribution to the sum in 
(1) of the amount 

■f{xu3 = X43). 



143a' 

After another "run" through the lattice sites, we will add a further contribution of 

Thus, continuing on this process ad infinitum, the overall contribution to the sum in (1) will 
simply be 



1 r 1 1 1 1 ^ 
a|i3 + T43 + 243 + 343 + ---P"^^^- 



It might seem that the coefficient of this term is simply the coefficient C43 that we seek to 
evaluate. However, we have a problem: this infinite sum is divergent, as can be seen from 
the fact that 1/43 > 1/100, 1/143 > 1/200, 1/243 > 1/300, and so on, so that 

1111 1111 



43 143 243 343 
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1 

^ 200 


+ 


1 


100 


300 


1 

100 




+ 


1 

3 + 



1 

4 + 



and this last sum is logarithmically divergent. 

What has gone wrong? Obviously, the correct coefficient cannot be infinite. The problem 
is that we have only considered one half of the derivative operation expressed in (1), namely, 
that for positive n. We have forgotten that, in travelling backwards along the lattice for 
negative n, and applying periodic boundary conditions at the boundary which we are taking 
to occur after n — —49, we will also "run past" the position at n = +43. Specifically, the 
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position n = —57 is, by the periodic boundary conditions, identically equivalent to n = +43, 
and this term in (1) will contribute an amount 



57a- 

Going back another N — 100 lattice sites, we will also find a contribution of 

and so on. Thus, we find that the true expression for C43 is simply 

_ If 1 1 1 1 1 1 -] 

~ a 143 ~57^143~157^243~257^'"J' ^ ' 

It is straightforward to see that this alternating series is now convergent; for example, by 
noting that 1/43 < 1/25, 1/57 < 1/50, 1/143 < 1/75, and so on, so that 

irl 1 1 1 1 If 111 lln2 

C43 < -\ \ h . . . > = \ 1 \ h...> = < 00. 

al25 50 75 100 / 25al 2 3 4 / 25a 

Now, you may have noted that, in "wrapping around" the N — 100 lattice sites in the 
above, the fact that N was even meant that the sign of the contribution to (1) in any "run" 
through the lattice (in one direction) was the same as any other "run" through the lattice 
(in that same direction). How would the situation change if were odd7 To examine this 
situation, let us consider the example of = 101, and let us continue to concentrate our 
attention on the position n — +43. Going through the positive-n terms in (1), we will pick 
up contributions of 



If 1 1 1 1 

al43~ 144^ 245 " 346 ^ 



■}, (7a) 



where the minus signs arise because the terms in (1) for even positive n come in with the 
minus signs. This sum is now, by itself, convergent. However, to obtain C43 we still need to 
add in the contributions from the ncgative-n terms in (1). The periodic boundary conditions 
now tell us that the position n = —58 is equivalent to n = 43 (because 43 + 58 = 101), and 
the coefficient of f{x — 58a) in (1) is positive, so we will obtain the contribution 

+ ■■■}, (7b) 



58 159 260 261 
which is also, by itself, convergent. 



4. An example of the second derivative 

The same process as outlined in the example above can be applied to the second-derivative 
operator (2) (and, indeed, any higher-order derivative that may be desired). Let us go back to 
the example oi N — 100, and consider again the position n — +43. Under periodic boundary 
conditions, the terms in (2) for positive n that contribution to this position are just 

^{4^ + Ti3^ + 2i3^ + 3^ + ■ • j^^^^^^- ^'^^ 
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Because these denominators are squared, rather than hnear, this sum converges on its own. 
The negative-n terms in (2) that contribute to this position are just 

+ 1^ + 2^ + 3^ + ■ ■ j^^^^^)' ^'^^ 

which also converges. For even N, these two sums contribute with the same sign. 

It is clear that the same of calculation would also apply for the case of n = +50, which 
is "half a lattice away" from n = 0, even though there is no corresponding n = —50 (which 
is, by the periodic boundary conditions, equivalent to n = +50); for this position, the two 
sums of the form (8a) and (8b) will be identical. 

Now, because the central {n = 0) term in the operator (2) is nonzero, and because each 
term for —n comes in with the same sign as that for +n, we also need to separately compute 
the coefficient cq under the assumption of periodic boundary conditions. For the case of 
N = 100, we simply obtain 

(9) 




1002 2002 3002 



where the factor of 2 in each term arises from the fact that we are adding in the contributions 
from both positive-n and negativc-n (which are both one lattice "revolution" away). 

Let us now go back to the odd-iV example of iV = 101. For the position n = +43, the 
positive-n terms that contribute are just 

2 f 1 1 1 1 



a2l432 1442 + 245^- 346^ + UN 
and the negative-n terms that contribute are just 

U'W^ + 1^ - 2^ + 3^ + ■ ■ ]^^'''^- 

The extra minus signs that arise from using an odd number of lattice sites make the overall 
sum converge more quickly than for an even value of N. For the "central" coefficient cq, we 
now obtain 

_2_r_7r_2 _2 2_ _2__ 1 

6^ 1012 2022 + 3032 ■■■/' ^ ^ 

where the minus signs are again due to the fact that we are alternately picking up the odd 
and even terms in (2). 



5. Exact, general expressions for the finite-lattice SLAC derivative operators 

Let us now make use of the examples provided in Sees. 3 and 4 to guide us in constructing, 
from first principles, the general expressions for the SLAC derivative operators on a finite 
lattice on which we assume periodic boundary conditions. 

Let us start with the first-derivative operator. We can write the expression (6) for C43 
with N = 100 in the form 

(lUoo ^ 1 y-f 1 ^ ] 

a tr'A 100k + 43 lOOA; + 57 / ' 
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where I am using the notation 



(m) iV 



for the coefficient of the ideal SLAC m-th derivative operator at position n for a finite lattice 
of N lattice sites. If we had chosen a general lattice position n (positive) rather than the 
particular choice of n — 43, we would have obtained 



(iLioo ^ \-^) 

n positive 

^ k=0 



I I 

frlllOO/c + n lOO/c + lOO-nJ' 



where the alternating sign out the front is due to the alternating sign in (1). Clearly, if we 
were to have a general (even) number of lattice sites, this result would become 



_ (-l)'*+i ^ 
n positive - - 2^ ] Nk + n ~ Nk + N - n 



(1) Neven 



k=0 

If we factorise out a factor of N from each of these denominators, we obtain 



n positive- k+l-n/Ny 



(1) N even 



Fortunately, the mathematicians have again evaluated this infinite sum for us, and it yields 
a remarkably simple result: 



k=0 



Thus, we find that 



(1) JVeven _ ( _-i)n+l cot( — ) 
C„posrt:ve-l iVa V iV j ' 



Let us now consider a negative value of n (but keeping N even). From (1), all of the terms 
come in with the opposite sign, but since (— 1)"""*"^ = (— 1)""*"^ and cot(— x) = — cot a;, we 
can see that the same expression can be used for negative n. Hence, in general (for even N), 
we have 

'■'<=S- = (-l)"«j^<=ot(^^). (13) 

It is clear that, for small n (namely, far from the application of the periodic boundary condi- 
tions), this formula approaches the infinite lattice result (1), because cotx = 1/ tanx ~ 1/x 
for small x, and hence 

(1) JVeven — ^ () 

On the other hand, near the edge of the application of periodic boundary conditions, namely, 
n ^ N/2, the absolute values of the coefficients (13) approach zero hnearly with {N/2 — n). 
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Let us now consider the case of an odd number of lattice sites A^. We again start at a 
position of positive n. The example (7a) shows us that the positive-n terms in (1) themselves 
contribute an amount that can be written 



-1 C M 1 2M 

t lini 2y — y 



101A; + 43 i' 



where we have written the two sums in this way (the first sum contributing to every other 
term in (7a), with the second sum contributing to every term) to maintain the closest possible 
relationship with the eyen-N analysis above. (The limiting process is necessary because each 
sum does not individually converge.) Likewise, by (7b) we see that the negative-n terms in 
(1) contribute an amount 



-1 r M 1 2M 

i lim 2V — r y.- 



101A; + 58 

Putting these together, we then have 

(iLioi ^ i lim i V ( = + = 1 - V 

^3 aM^oo\^QV202A; + 43 202A; + 58; 



+ 



J01A; + 43 101A; + 58. 
Generalising to the case of general (odd) A^, and positive n, we clearly obtain 

^c^^... = t^hm ^.^.^ + 



-t\n+l ( ^ / 

lim V( 



■n positive ^ M^oo\{^\2Nk + n 2Nk + N -n 



2M 



fc=0 

Again factorising A^ from these denominators, we then obtain 



(_l)n+i ( M 
n positive = M-^^[^^\k + n/2N ' k + l/2-n/2N^ 



(l)^Arodd 



lim ^ V — - + 



2Af / 2 I 



k=0 



k + n/N k + l-n/N^ 



Fortunately, the mathematicians have also evaluated this combination of infinite limits of 
sums for us, and again it is a remarkably simple result: 



M^^X + n/2N + A; + 1/2 - n/2 AT ) ' 



1 

+ 



/ 7rn\ 



i/2N k+l/2-n/2N ) ^^\k + n/N k + l-n/N)] (14) 

= TT CSC I 

We thus find that 

positive — I ^) '^i>(.\^ ^ 



As before, this expression is odd under n — > —n, in accordance with (1), and so for a general 
position n (for odd N) we have the result 

'"ir = (-l)""l^csc('^). (15) 



Again, since esc a; = 1/ sin a; ^ 1/x for small x, we find that far from the application of 
periodic boundary conditions we regain the infinite- lattice result: 

(i)jvodd {-iT^^ — ^0 

""^^^ na N 

However, in contrast to the case of even iV, for odd N we find that, near the edge of the 
application of periodic boundary conditions, namely, n N/2, the absolute value of ^^^c^^^'^ 
approaches the finite value n/Na. If it weren't for the periodic boundary conditions, the 
absolute value of the coefficient function at the limit of this boundary (between lattice sites) 
would only have been 1/na = 2/Na; thus, the applications of periodic boundary conditions 
for odd N has actually increased the magnitude of the coefficient function at the boundary 
by a factor of n/2, or around a 57% increase. Indeed, this is a general phenomenon: 



(1) JVodd 



1 

> — , 

na 



whereas 

1 



{l)„Neven 



< . 

na ' 



although in each case the magnitude of the inequality is relatively mild, as we shall see 
shortly in Sec. 6. 

Let us now turn to the second-derivative operator. The expressions (8a) and (8b) clearly 
combine to give us 



r\ OO 
(2)„100 _ _f_ 



-43 



a~oi(100A; + 43)2 {lOOk + 57y 



Going to general (even) and general n (positive or negative, since they come in with 
the same sign for the second-derivative), we clearly have 



(2) N even 



^1 (' ATI- _l_ ^^2 



a2 ^^[{Nk + n)^ ' {Nk + N-nfl 
and factorising out N"^ from the denominators we have 

(2) Afeven _ 2( — 1)"+^ ^ j ^ . \ 

'^^^ {Naf {k + n/Nf {k + I - n/Nf 

Again, the mathematicians give us a magic result for the infinite sum: 
" ' 1 1 1 „ jTm\ 



£1 (^W^ + (^ + 1 - n/Nf] = AT j ' ^1^) 
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so that 

(^)c-- = 2(-l)'^+^(^ Jcsc2(^) . (17) 

From the way we have constructed it, this expression clearly also applies to the (single) 
position n = +N/2, that is half a lattice away from the position n = at which we are 
applying the derivative. Clearly, far from the application of periodic boundary conditions, 
we again regain from (17) the infinite-lattice result (2). 
For n — 0, the expression (9) can be written 



(2)JVeven _ _'^\'^ _L V ± I 



and since 



°° 1 7r2 



we find 



(2)JVeven _ , ^\ Qg^ 

It can be shown that the additional contribution here for finite N is equal to the difference 
between the finite-lattice result (17) and the infinite-lattice result (2) forn ^ if we assume 
the expression to hold true for real n and take the limit n ^ 0. 

Finally, we consider the case of odd N. The sum in (10a) can be written 



2 



aH ,t^o (202A; + 43)2 (lOl/c + 43)^ j ' 

where in this case we don't need to be careful about the limiting process because each sum 
is separately convergent. Likewise, the sum in (10b) can be written 

2 r oo 2 °° 1 ] 

~a2 f ,5 (202A; + 58)2 " E (lou + 58)2 |' 

Combining these expressions, we then have 

(2) 101 ^ ^ J 2 2 1 1 



(202A; + 43)2 (202A; + 58)2 (101A; + 43)2 (lOlA; + 58)2 j ' 

For general (odd) and n 7^ 0, this becomes 

00 

'C'^Jn" = ^ '- ^< 

1 



(2)^Arodd ^ 2(-l)"+^ j 2 2 



{Nk + n)2 {Nk + N- nf 
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and on factorising out N"^ we get 



(2LiVodd ^ 2(-ir+i - f 1/2 1/2 

[Naf ^^\{k + n/2NY {k + l/2-n/2NY 

1 1 

+ 



{k + n/Nf {k + l-n/Nf 
Again, the mathematicians have evaluated this sum for us: 
1/2 1/2 



k=0 



SO that 



{k + n/2Ny {k + l/2-n/2Nf 



1 1 1 „ /7rn\ /7rn\ 

+ 77^-; TT^^ = ^ CSC — cot — 1 



{k + n/Nf {k + l-n/Nf \ \NJ \NJ' 



(2) JVodd 



For n = 0, the expression (11) can be written 




and so, using (5), we find 



(2) iVodd _ _]^( . _ J_\ 



The expressions (20) and (21) are in agreement with those hsted exphcitly by DreU, Weinstein 
and Yankielowicz [2] in the first of their seminal papers. Note, however, that these authors 
only considered the case of an odd number of lattice sites (which is written as 2A'" + 1 in their 
notation); and as no explicit expressions for the first-derivative operation in position space 
are provided in [2], further comparisons with the expressions listed above are not possible. 

Clearly, the expressions above for the second-derivative operator on finite lattices also 
reduce to the infinite-lattice results for n/A^ — > 0, as was the case for the first-derivative 
operators. Near the application of periodic boundary conditions, the even-A" expression 
(17) approaches a finite limit, whereas the odd-A^ expression (20) approaches zero, which 
is a reversal of roles from the case of the first-derivative operators (for which the even-A" 
coefficient vanished linearly, and the odd-A^ coefficient approached a finite value). 



6. Stochastic implementation of the finite-lattice SLAC derivative operators 

Let us now turn to the reason why we are considering the practical implementation of 
the SLAC derivative operators at all, namely, my proposal in [1] that they be implemented 
in a stochastic fashion. How must this proposal be modified to accommodate the corrections 
obtained in Sec. 5 to incorporate exactly the implementation of periodic boundary conditions 
for a finite lattice? 
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I propose the following: that the probabilistic weight assigned to the computation of any 
difference (sum) in the first (second) SLAC derivative operator at distance n be calculated 
using the same value 1/n (2/n^) as was proposed in [1] on the basis of the truncation of 
the infinite-lattice operators, but that the result of such sum (difference) — if calculated — be 
multiplied by a "correction factor" that incorporates exactly the effects of restricting the 
formalism to a finite lattice and applying periodic boundary conditions. 

My reasons for proposing this solution are twofold. Firstly, as noted in Sec. 5 (and as will 
become explicit below), the required "correction factors" are, in all cases, simply numbers 
of order unity, which do not affect the fluctuation properties of the proposed stochastic 
implementation in any functionally critical fashion; we simply weight some positions a little 
more, or less, than for the case of an infinite lattice. Secondly, it is clear that if this is 
all that we need to do, then it would be silly to tamper with the extraordinarily efficient 
method of random selection described in [1] for probabihstic weight functions of 1/n or 
because the decision to compute or not compute needs to be made for every distance n for 
every application of the derivative operator, and may, if blown out to a set of floating-point 
computations, become the time-critical part of the algorithm that would limit its practical 
implementation on a digital computer. (It is important to note that the floating-point 
multiplication of the result of the sum or difference by the correction factor need only be 
done if the sum or difference is actually computed, and hence the average number of extra 
floating-point operations according to this proposed solution is itself proportional to InN, 
rather than being proportional to if the probabilistic weights were instead modified.) 

Prom the results obtained in Sec. 5, it is simple to write down expressions for these 
correction factors. If we define 

(m) JV 

to be the correction factor for the distance-n sum or difference for the m-th order SLAC 
derivative operator on a (one-dimensional) lattice of N sites, then clearly 




From the results (13), (15), (17) and (20) we then simply have 




(22a) 
(22b) 
(22c) 
(22d) 



where 



A^ 



Tin 




(23) 



and from (18) and (21) we have 




2 



(24a) 



(2) iVodd 




(24b) 
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As noted above, the expressions (22) are relatively "boring" functions of a. They all approach 
unity for a — > 0. For a. — > 7r/2, the expressions (22a) and (22d), involving cot a, fall smoothly 
to zero, whereas the expressions (22b) and (22c) rise to the finite values 7r/2 and 77^4 
respectively. 

7. Conclusions 

We have found that it is relatively easy to get a good physical understanding of the SLAC 

derivative operators on finite one-dimensional lattices, in position space, by considering 
explicitly the imposition of periodic boundary conditions on the corresponding expressions 
for infinite lattices. In all cases this leads to small adjustments of coefficients, which can 
be written down in a very simple form. There are qualitative differences between these 
correction factors depending on whether the lattice has an odd or an even number of sites, 
but in all cases the factors arc numbers of order unity. This means that the stochastic 
implementation of the operators proposed in [1] requires only superficial changes in order to 
make the operators completely optimal for any finite lattice. 
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